ADP-02-74-T513 
JLAB-THY-02-26 

Moments of Isovector Quark Distributions from Lattice QCD 

W. Detmold, 1 W. Melnitchouk, 2 and A. W. Thomas 1 

1 Special Research Centre for the Subatomic Structure of 
Matter and Department of Physics and Mathematical Physics, 
University of Adelaide, SA 5005, Australia 
2 Jefferson Lab, 12000 Jefferson Avenue, Newport News, VA 23606 

g i Abstract 

\ We present a complete analysis of the chiral extrapolation of lattice moments of all twist-2 



isovector quark distributions, including corrections from Nir and An loops. Even though the A 
resonance formally gives rise to higher order non-analytic structure, the coefficients of the higher 
order terms for the helicity and transversity moments are large and cancel much of the curvature 
| generated by the wave function renormalization. The net effect is that, whereas the unpolarized 
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<Q ■ moments exhibit considerable curvature, the polarized moments show little deviation from linearity 

. as the chiral limit is approached. 
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I. INTRODUCTION 

Resolving the quark and gluon structure of the nucleon remains one of the central chal- 
lenges in strong interaction physics JT|. Information about the nucleon's internal structure 
is parameterized in the form of leading twist parton distribution functions (PDFs), which 
are interpreted as probability distributions for finding specific partons (quarks, antiquarks, 
gluons) in the nucleon in the infinite momentum frame. PDFs have been measured in a 
variety of high energy processes ranging from deep-inelastic lepton scattering to Drell-Yan 
and massive vector boson production in hadron-hadron collisions. A wealth of experimental 
information now exists on spin-averaged PDFs, and an increasing amount of data is being 
accumulated on spin- dependent PDFs @. 

The fact that such a vast array of high energy data can be analyzed in terms of a universal 
set of PDFs stems from the factorization property of high energy scattering processes, in 
which the short and long distance components of scattering amplitudes can be separated ac- 
cording to a well-defined procedure. Factorization theorems allow a given differential cross 
section, or structure function, F, to be written (as a function of the light-cone momen- 
tum fraction x at a scale Q 2 ) in terms of a convolution of hard, perturbatively calculable 
coefficient functions, Cj, with the PDFs, fi, describing the soft, non-perturbative physics 



where /x is the factorization scale. The coefficient functions are scale and process dependent, 
while the PDFs are process independent, and hence can be used to parameterize a wide 
variety of high energy data. 

Because the PDFs cannot be calculated within perturbative QCD, the approach com- 
monly used in global analyses of high energy data is to simply parameterize the PDFs, 
without attempting to assess their dynamical origin ||, [], Once fitted at a particular 
scale, they can be evolved to any other scale through the DGLAP Q 2 -evolution equations |§. 
The focus in this approach is not so much on understanding the non-perturbative (confine- 
ment) physics responsible for the specific features of the PDFs, but rather on understanding 
the higher order QCD corrections for high energy processes. 

In a more ambitious approach one would like to extract information about non- 
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perturbative hadron structure from the PDFs. However, without an analytic solution of 
QCD in the low energy realm one must rely to varying degrees on models of (or approx- 
imations to) QCD within which to interpret the data. An extensive phenomenology has 
been developed over the years within studies of QCD-motivated models, and in some cases 
remarkable predictions have been made from the insight gained into the non-perturbative 
structure of the nucleon. An example is the d—u asymmetry, predicted || on the basis of the 
nucleon's pion cloud Jl^|, which has been spectacularly confirmed in recent experiments at 
CERN, Fermilab and DESY f I T]| . Other predictions, such as asymmetries between strange 
and anti-strange [|l^], and spin-dependent sea quark distributions, Au — Ad Hl3f , still await 
definitive experimental confirmation. Note that none of these could be anticipated without 
insight into the non-perturbative structure of QCD. 

Despite the phenomenological successes in correlating deep-inelastic and other high energy 
data with low energy hadron structure, the ad hoc nature of some of the assumptions made 
in deriving the low energy models from QCD leaves open questions about the ability to 
reliably assign systematic errors to the model predictions. One approach in which structure 
functions can be studied systematically from first principles, and which at the same time 
allows one to search for and identify the relevant low energy QCD degrees of freedom, is 
lattice QCD. 

Lattice QCD is rapidly developing into an extremely useful and practical tool with which 
to study hadronic structure |T3. There, the equations of motion are solved numerically 
by discretizing space-time into a lattice, with quarks occupying the lattice sites and gluons 
represented by the links between the sites. Meaningful numerical results can be obtained 
by Wick rotating the QCD action into Euclidean space. However, because the leading twist 
PDFs are light-cone correlation functions (involving currents with space-time separation, 
z 2 — (ct) 2 ~ 0), it is, in practice, not possible to calculate PDFs directly in Euclidean space 
- in this case a null vector would require each space-time component to approach zero. 
Instead, one uses the operator product expansion to formally express the moments of the 
PDFs in terms of hadronic matrix elements of local operators, which can then be calculated 
numerically. 

In relating the lattice moments to experiment, a number of extrapolations must be per- 
formed. Since lattice calculations are performed at some finite lattice spacing, a, the results 
must be extrapolated to the continuum limit, a — > 0, which can be done by calculating 
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at two or more values of a. Furthermore, finite volume effects associated with the size of 
the lattice must be controlled — working with a volume that is too small can result in 
the omission of important physics, arising from the long-range part of the nucleon wave 
function. Finally, since current lattice simulations typically use quark masses m' att above 
30 MeV, an extrapolation to physical masses, mj hys « 5 MeV, is necessary. Earlier work on 
moments of spin-averaged PDFs [ 15 , P^[ ] found that whereas the lattice calculations yielded 
results typically 50% larger than experiment when extrapolated linearly to mj hys , inclusion 
of the nonlinear, non-analytic dependence on m q arising from the long range structure of 
the nucleon removes most of the discrepancy. 

In this paper we extend the analysis to the polarized sector, which is important for 
several reasons. Firstly, for many years lattice calculations of the axial vector charge, g^, 
have tended to lie 10% or more below the experimental value determined from neutron (3 
decay. Since this represents one of the benchmark calculations in lattice QCD, it is vital 
that the source of this discrepancy be identified. A preliminary analysis of the effects of 
chiral loops found [17], pl| that the inclusion of the leading non-analytic (LNA) behavior 
associated with nN intermediate states in the extrapolation of qa to mj hys decreased the 
value of gAi thereby making the disagreement worse. On the other hand, one knows that the 
A resonance plays an important role in hadronic physics, so a more thorough investigation of 
its effects on spin-dependent PDFs is necessary before definitive conclusions can be drawn. 
Indeed, we find that although the A contributions formally enter at higher order in m n , the 
coefficients of these next-to-leading non-analytic terms are large, and their effects cannot be 
ignored in any quantitative analysis. In addition, since there are currently no data at all on 
the transversity distribution in the nucleon, lattice calculations of several low transversity 
moments provide predictions which can be tested by future measurements. In order to make 
these predictions reliable, it is essential that the lattice calculations be reanalyzed to take 
into account the chiral corrections entering extrapolations to m^ hys . 

The remainder of this manuscript is structured in the following manner. In Section || we 
describe the calculations of the moments of PDFs from matrix elements of local operators, 
and summarize the details of extant lattice calculations. In Section |TTTj we first examine the 
constraints from chiral perturbation theory and the heavy quark limit on the behavior of 
the moments of the various distributions as a function of the quark mass. The importance 
of higher order terms in the chiral expansion is then investigated within a model which 
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preserves the non-analytic behavior of chiral perturbation theory. This information is used 
to construct effective parameterizations of the quark mass dependence of these moments, 
which are then used to extrapolate the available lattice data in Section |TV| . Finally, in 
Section [V] we discuss the results of this analysis and draw conclusions. 

II. LATTICE MOMENTS OF PARTON DISTRIBUTION FUNCTIONS 
A. Definitions 

The moments of the spin-independent, q = q^+q, helicity, Aq = q^ — q^, and transversity, 
5q = q T — q ± , distributions are defined as: 

(x n ) g = [ dxx n [q(x) - {-l) n q{x)} , (2a) 
Jo 

(x n ) Aq = [ dxx n [Aq(x) + (-l) n Aq(x)} , (2b) 
Jo 

(x n ) Sq = [ dxx n [5q(x)-{-l) n 5q(x)), (2c) 
Jo 

where g^U) corresponds to quarks with helicity aligned (anti-aligned) with that of a longi- 
tudinally polarized target, and g 1 "^ corresponds to quarks with spin aligned (anti-aligned) 

with that of transversely polarized target 1 . At leading twist, these moments depend on 



ground state nucleon matrix elements of the operators 

Qi^-ixn = { n-\ ^ 7 {mi £>M2 . . . D ^}^ ; ( 3a ) 

= i n ~ x i> 7 {M1 75 ■ ■ ■ D"»ty , (3b) 

Og* 1 "*"" = i n ~ x i> <r a{ ^ 75 P M2 ■ ■ ■ D^ty , (3c) 

respectively. Thus, for a nucleon of mass M, momentum P, and spin S, one has: 

(P,S\0^-^\P,S) = 2 pte ...pM - traces, (4a) 

(P, S|C^- M "|P, S) = 2 (x^Aq M S { ^P^ ■ ■ ■ P^" } - traces , (4b) 

(P,S\O a & ^ n \P,S) = 2 (x n -% M S [a P^ l] P^ ---P^ - traces, (4c) 



Note that from their definition, Eqs. (g), the moments alternate between the total (q + q) and valence 
(q — q) distributions, depending on whether n is even or odd. 
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where the braces, { • • • }([■■•]), imply symmetrization (anti-symmetrization) of indices, 
and the 'traces' (containing contractions <7 WW , etc.) are subtracted to make the matrix 
elements traceless in order that they transform irreducibly. At higher twist (suppressed 
by powers of 1/Q 2 ), more complicated operators involving both quark and gluon fields 
contribute. 

B. Lattice Operators 

The construction of the relations between moments of PDFs and matrix elements of 
local operators relies on the symmetry group of the Euclidean space in which one works. 
When formulated on a discrete space-time lattice, the symmetry group is reduced and the 
discretized implementation of these operators introduces several technical complications. 

The discrete nature of the lattice topology means that the symmetry group of the Eu- 
clidean continuum, the orthogonal group 0(4), is broken to its hyper-cubic subgroup, H(4) 
(the group of 192 discrete rotations which map the lattice onto itself) Unfortunately, 
operators in irreducible representations of 0(4) may transform reducibly under H(4) and 
this may result in mixing with operators from lower dimensional multiplets under renor- 
malization. Consequently, care must be exercised in the choice of operators on the lattice. 
For example, the continuum operator 0%", which corresponds to the momentum carried by 
quarks, can be represented on the lattice by either O^f 1 = ipj^D^i/j (belonging to a 6 repre- 
sentation) or Oq = ifj^D^ip — | X/i=i ^7*-^V (belonging to a 3 representation). This may 
be regarded as an advantage since in the a — > limit these operators are identical and any 
difference at non-zero lattice spacing allows for an estimate of the remaining finite size lattice 
artifacts to be made. In practice, this is currently somewhat ambitious, as the operator O g a ^ 
requires that the hadron source should have non-zero momentum components, which leads 
to a statistically less well determined result. Consequently, for the operator 0^ u we retain 
only the data corresponding to O q . The operators associated with the unpolarized n = 2 
and 3 moments are given by 0\ llA} - \ £t 2 Of 4} and Of 44} + 33} - O g {1133} - O g {2244} , 
respectively. 

For the spin-dependent moments, the operator corresponding to the axial charge is given 
by 0\ q = ip^j 3 ^. However, for the n = 1 moment one can have on the lattice either 
£>2 = ^ 7 5 7 { 1 D 3 >V' or = ^7V 3 £> 4} ^- Once again, since the operator 0£ requires 
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non-zero momentum, we shall keep only the data corresponding to the better determined 
O^j operator. The operators required to calculate other moments in Eqs. (0) are described 



in Ref . 20 



For spin greater than 3, there are no unique, irreducible representations in H(4) for the 
twist-2 operators. This means that the operators for moments n > 3 will inevitably mix with 
lower dimensional (or lower spin) operators. To unambiguously extract information about 
these moments, one would need to consider all representations for a given spin, and, with 
sufficiently accurate data, deduce the matrix elements of the high spin operators from the 
low spin operators with which they mix. Because of these difficulties, all lattice calculations 
have so far been restricted to moments with n < 3. Nevertheless, some features of the PDFs 



can be reasonably reconstructed from just the lowest few moments, as described in Ref. [16 



Further subtleties arise when we consider the non-perturbative renormalization of these 
operators and their matching to other renormalization schemes. An operator, Oi att (a), 
calculated using the lattice regularization scheme, is connected to other schemes, for example 
MS, by a renormalization factor: 

CWM = Z (jJi,a) Ciatt(a) , 

where /x is the renormalization scale. To provide results in standard schemes, the renormal- 
ization functions, Z@, must therefore be calculated for each operator used. While this is 
done perturbatively in most calculations, non-perturbative determinations now exist |2T| . In 
what follows, results are presented in the MS scheme at a renormalization scale /i 2 ~ 4 GeV 2 . 



Lattice Calculations 



The first calculations of structure functions within lattice QCD were performed in the 
late 1980s by Martinelli and Sachrajda. Their pioneering calculations of quark distributions 
of the pion []22| and nucleon [^3| were ambitious, given the speed of the computers available 
at the time. More recently, various calculations of greater precision have been performed 

j, [34], [35|, 36, 37]. In the present analysis we will 



191, 20, 24, 25, 26, 27, G8L E9l ROL Bl 



focus mainly on the more recent QCDSF |30], |3I 



331 and LHPC-SESAM simulations 19 



20| , p4| . The older data sets from Gupta et al. [[35] have large uncertainties associated with 
renormalization, while the statistical precision of Refs. [25], is comparatively low. In 
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addition, several groups (notably the KEK |7|, Riken-BNL-Columbia (RBC) |§, |37 



and SESAM |29|] collaborations) have put particular emphasis on the n = moments of 
the helicity and transversity distributions — the axial and tensor charges. The simulations 
have been made using various quark and gluon actions, on different lattices and at different 
couplings. They have been performed primarily in the quenched approximation, although 
more recently the LHPC |2IJ and UKQCD/QCDSF groups have begun to investigate 
the effects of unquenching. In Table | we summarize the data used here, for which PDF 
moments and the corresponding pion masses are published. 

Before including the data sets in our analysis, we impose a simple cut to reduce finite 
volume effects. In lattice calculations of any observable, one must ensure that the lattice size 
is large enough for results not to be dependent on the (unphysical) boundary conditions. 
This applies particularly to calculations involving low energy states such as the nucleon 
where the effects of the pion cloud are known to be especially important. Being the lightest, 
and therefore longest range, asymptotic correlation of quarks and gluons, pions are most 
sensitive to the boundary conditions. To avoid these difficulties, we require that the lattice 
volume is large enough that a pion will fit comfortably within it without "feeling the edges of 
the box". A pion of mass m n has a corresponding Compton wavelength of order ~ l/m^, 
and, to avoid interference of the pion with its periodic copies, we require that the smallest 
dimension of the lattice box (L) satisfies the constraint 

L > A\ n ~ — . (5) 



The factor of 4 in this formula is popular fL9 |, although somewhat arbitrary. This argument 
indicates that the lowest mass data point of Ref. |H[ and the lightest unquenched points 
from Ref. |2(J should be excluded from the analysis. 

In terms of quark flow, the evaluation of matrix elements of the operators in Eqs. @ 
includes both connected and disconnected diagrams, corresponding to operator insertions 
in quark lines which are connected or disconnected (except through gluon lines) with the 
nucleon source — see Fig. |I|. Since the numerical computation of disconnected diagrams is 
considerably more difficult, only exploratory studies of these have thus far been completed 
p9| , and the data analyzed here include only connected contributions. However, because 
the disconnected contributions are flavor independent (for equal u and d quark masses), 
they exactly cancel in the difference of u and d moments. Therefore, until more complete 
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(a) (b) 

FIG. 1: Connected (a) and disconnected (b) contributions to the matrix elements of an operator 
(indicated by the cross). Such diagrams occur in quenched QCD as well as in full QCD. 

lattice simulations become available, one can only compare lattice moments of the flavor 
non-singlet u — d distribution with moments of phenomenological PDFs H]. 

Whilst the chiral behavior of quenched QCD is different from that of full QCD P5[ , 
in the region where current lattice data from both quenched and unquenched simulations 
are available, the differences are well within the statistical errors, indicating that internal 
quark loops do not play a significant role over this mass range. As calculations begin to 
probe lighter quark masses, the differences should become more apparent and it will become 
necessary to analyze quenched and unquenched data separately |]39" |. Until the differences 



become statistically distinguishable, however, we shall combine the data from the two sets 
of simulations. 



III. CHIRAL BEHAVIOR OF PDF MOMENTS 

To compare the lattice results with the experimentally measured moments, one must 
extrapolate the data from the lowest quark mass used (~ 50 MeV) to the physical value 
(~ 5 MeV). This is commonly done by assuming that the moments depend linearly on the 
quark mass. However, as discussed in Ref. [0, such a linear extrapolation overestimates 
the experimental values of the unpolarized isovector moments ||, [|, |6| by some 50% in all 



cases. Since the discrepancy persists in unquenched simulations [0, |20], gj] , it suggests that 
important physics is being omitted from either the lattice calculations or their extrapolations. 
In Refs. |15|, the chiral behavior of the moments of the unpolarized isovector distributions 
was found to be vital in resolving this discrepancy. Here we summarize the results of the 
earlier unpolarized study, and extend the analysis to the moments of the spin-dependent 
isovector distributions in the nucleon. 
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Reference 


Q/U 


Quark Action 


Lattice 


a [fm] 


[GeV] 


Moments 


Symbol 


QCDSF [30] 


Q 


Wilson 


16 3 x 32 


0.1 


0.6 - 1.0 


All 


▲ 


QCDSF [31] 


Q 


Wilson 


24 3 x 32 


0.1 


0.35 - 0.6 


All 


■ 


QCDSF [32] 

v L J 


Q 


NPIC 


16 3 x 32 


0.1 


0.6 - 1.0 


All 


x 


QCDSF [33] 


Q 


NPIC 


16 3 x 32 


0.1 


0.65 - 1.2 


(l>Ag, (x 2 )Aq 


♦ 




Q 


NPIC 


24 3 x 48 


0.075 


0.7 - 1.2 


(1)A ?I {x 2 )Ag 


• 




Q 


NPIC 


32 3 x 48 


0.05 


0.6 - 1.25 


(1)A ?I {x 2 )Ag 




MIT [20] 


Q 


Wilson 


16 3 x 32 


0.1 


0.58 - 0.82 


All 


★ 


MIT-SESAM [20] 


u 


Wilson 


16 3 x 32 


0.1 


0.63 - 1.0 


All 





MIT-SCRI [20] 


u 


Wilson 


16 3 x 32 


0.1 


0.48 - 0.67 


All 


□ 


D^EK [27, 2SL |_. 


□ & 


Wilson 


16 3 x 20 


0.14 


0.52 - 0.97 


(l>Ag, (l)<yg 


▼ 



TABLE I: Simulation parameters for lattice calculations of the moments of PDFs included in our analysis. Q/U corresponds to 
quenched/unquenched simulations, and NPIC denotes the nonperturbatively improved clover quark action. "All" moments correspond 
to (x l ) q for i = 1,2, 3, {x l )& q for i = 0, 1, 2, and {x l )s q for i = 0,1. The symbols shown in the final column correspond to those plotted in 
Figs. 8, 9 and 10. 



A. Chiral Symmetry and Leading Non-analytic Behavior 



The spontaneous breaking of the chiral SUl(N/)xSUr(N/) symmetry of QCD generates 
the nearly massless Goldstone bosons (pions), whose importance in hadron structure is well 
documented. At small pion masses, hadronic observables can be systematically expanded 
in a series in m n — chiral perturbation theory (%PT) [|4(|. The expansion coefficients are 
generally free parameters, determined from phenomenology. One of the unique consequences 
of pion loops, however, is the appearance of non-analytic behavior in the quark mass. From 
the Gell-Mann-Oakes-Renner relation one finds that m 2 ~ m q at small m^, so that terms 
involving odd powers or logarithms of m n are non-analytic functions of the quark mass. 
Their presence can lead to highly nonlinear behavior near the chiral limit (m n — > 0) . 



Because the non-analytic terms arise from the infrared behavior of the chiral loops, they are 
generally model independent. 

The leading order (in m n ) non-analytic term in the expansion of the moments of PDFs 
was shown in Ref . |42|] to have the generic behavior m 2 log arising from 7r N intermediate 
states. This was later confirmed in HU, where the coefficients of these terms were 



also calculated. In Ref. a low order chiral expansion for the moments of the non-singlet 
distribution, u — d, was developed, incorporating the LNA behavior of the moments as a 
function of m q and also connecting to the heavy quark limit (in which quark distributions 
become (^-functions centered at x — 1/3) [jnj. For the moments of the unpolarized isovector 



distribution, these considerations lead to the following functional form for the moments [16 

1 + c LNA ml log 77 + b n - , (6) 
ml + fi 2 J ml + m 2 bn 

where (for n > 0) the chiral coefficient clna = —(1 + 3g\)/ (4:irf n ) 2 [f|3|], and b n is a constant 
constrained by the heavy quark limit: 

K = ^ - a n (l - /i 2 c L NA) • (7) 

The n = moment, which corresponds to isospin charge, is not renormalized by pion loops. 
The parameter /i is introduced to suppress the rapid variation of the logarithm for pion 
masses away from the chiral limit where %PT breaks down. Physically it is related to the 
size of the nucleon core, which acts as the source of the pion field ||41|| . Finally, the fits to 
the data are quite insensitive to the choice of m bn (as long as it is large), and it has been 
set to 5 GeV for all n gTf . 
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A similar analysis leads to analogous lowest order LNA parameterizations of the mass 



dependence of the spin-dependent moments [IS 



2 \ 2 

mi \ „ , mi 



(x n } A u-Ad = Ao n 1 + Ac L NA?r4log f + Ab n n , (8) 

and 

(Tfl^ \ 771 ^ 

1 + Sc LNA ml log ^ + / , (9) 

where the LNA coefficients are given by Aclna = —(1 + 2g^)/(47r/ 7r ) 2 and <5clna — — (1 + 
4g\) I '[2(47r A) 2 ] O- In the heavy quark limit, both Au(x) — Ad(x) and Su(x) — 5d(x) are 
given by — 1/3) [Q, which leads to the constraints: 

A&n = - Aa n (1 - /i 2 Ac LNA ) , (10) 

and 

5 

5b n = - 5a n (l- ffdcLNA) ■ (11) 

These are the most general lowest order parameterizations of the twist-2 PDF moments 
consistent with chiral symmetry and the heavy quark limits of QCD. 



B. Phenomenological Constraints 



In Refs. [|T^, [L6| we presented analyses of unpolarized data based on Eq. @, where it 
was concluded that current lattice data alone do not sufficiently constrain the extrapolation 
of these moments, and more accurate data at smaller quark masses (m q < 15 — 20 MeV) 
are required to determine the parameter /i. In that work, a central value of /x = 500 MeV 
(550 MeV when the heavy quark limit was not included |[15|| ) was chosen as it best reproduced 
both the lattice data and the phenomenological values at the physical point. However, the 
systematic error on this parameter is very large; indeed, the raw lattice data are consistent 
with /i = (a linear extrapolation). 

In order to make the phenomenological constraint of \i more quantitative, we employ the 
following measure of the goodness of fit of the extrapolated values (at mf ys ) of the first three 
non-trivial unpolarized moments to the phenomenological values, (x*)^^, as a function of 
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FIG. 2: The goodness of fit of the extrapolated values of the first three non-trivial moments to the 



phenomenological values as a function of \i calculated using Eq. (12). 



//: 



, CxptN 2 



(12) 



We assume that both the lattice data for the unpolarized moments and their extrapolation 
based on Eq. (||) are correct, and use the phenomenological values for these moments as a 
constraint. 

The behavior of the function xil 1 ) is shown in Fig. and the best value of \i is indeed 
found to be 500 MeV. This value is also comparable to the scale at which the behavior 
found in other observables, such as magnetic moments and masses, switches from smooth 
and constituent quark-like (slowly varying with respect to the current quark mass) to rapidly 
varying and dominated by Goldstone boson loops. For fits to lattice data on hadron masses, 



Leinweber et al. found values in the range 450 to 660 MeV |45fl when a sharp momentum cut- 
off was used. The similarity of these scales for the various observables is not coincidental, but 
simply reflects the common scale at which the Compton wavelength of the pion becomes 
comparable to the size of the bare nucleon. The value of /i is also similar to the scale 
predicted by the x 2 fits to the model discussed in the following section (see also Ref. HI 
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C. A Intermediate States 



When we come to the calculation of polarized PDFs, there is considerable phenomeno- 
logical evidence to suggest that the A resonance will play an important role. Within the 
cloudy bag model (CBM), a convergent perturbative expansion of the physical nucleon, in 



terms of the number of virtual pions, required the explicit inclusion of the A-isobar [BS, |47 



- see also Ref. j|8| for a recent, fully relativistic investigation. Without the A, the ratio of 
the bare to renormalized pion-nucleon coupling constant was found to be very large (as in 
the old Chew- Wick model). With it, they typically agree to within 10-15%. The essential 
physics is that the vertex renormalization associated with coupling to the A or to an N-A 
transition compensates almost exactly for the reduction caused by wave function renormal- 
ization. Of course, the same mechanisms apply to the renormalization of the axial charge, 
gAi as to the pion nucleon coupling, g n NN- 

In the limit that the A is degenerate with the nucleon, AM = — M — > 0, the leading 
non-analytic contribution from the A is of the same order as that arising from the nucleon, 
namely m^logm^. In the limit AM — > oo, the A contribution can be integrated out, and 
it formally does not make any non-analytic contribution. For a finite, but non-zero AM, 
the vertex renormalization involving the A is not a leading non-analytic term, but instead 
enters as m^j AM. However, the coefficient of this next-to-leading non-analytic (NLNA) 
term is huge |49| - roughly three times bigger than the m\ term in the expansion of the 



nucleon mass Faced with such a large coefficient, one cannot rely on naive ordering 

arguments alone to identify the important physics. 

The solution adopted by Leinweber et al. f|5] in the analysis of the chiral behavior of 
baryon masses was to calculate corrections arising from those pion loop diagrams respon- 
sible for the most rapid variation with m q . The finite spatial extension of the pion source 
leads naturally to an ultraviolet cut-off at the 7rNN and 7riVA vertices ||50|| . The parame- 
ter, A ~ 1/R (with R the size of the source) associated with these vertices is constrained 
phenomenologically. This approach ensures that the LNA and NLNA behavior of is 
reproduced in the — > limit, while the transition to the heavy quark limit (m, > A), 
where pion loops are suppressed as inverse powers of m^, is also guaranteed. Alternatively, 
one can study the variation of PDF moments with m n within a model, such as the cloudy 
bag H7], |5i"|], which also ensures the full LNA and NLNA behavior of x?T, and in addition 
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provides a simple physical interpretation of the short-distance contributions (in this case 
through the MIT bag model). Rather than rely on a specific model, in the present analy- 
sis we adopt the approach of Ref. and calculate the pion loop integrals with hadronic 
vertices constrained phenomenologically. 

The overall renormalization of the forward matrix elements of the operators of Eqs. ([3]) 
in nucleon states is then given by: 

(iV|Of -^|iV) dressed = ^(N\Of" n \N) haxe , i = q,Aq,5q , (13) 

where Z 2 is the wave function renormalization constant, 

Z-^l + Zf + Z 2 A , (14) 

and Zi are the vertex renormalization constants described below. The iV and A contributions 
to the wave function renormalization, illustrated in the first row of Fig. |||, are given in the 
heavy baryon limit 2 by [[5l|j : 

rN 39 2 a r k A u\k)dk 



2 " WKyJo ~^T' (15) 

7 a 4 ^ A 3g\ [<*> fcV(fc)dfc 
2 *& NN (±KU) 2 h ^)(^) + AM) 2 ' 1 ) 



where u(k) = \Jk 2 + m 2 is the pion energy, and u(k) is the form factor parameterizing the 
momentum dependence of the ttNN and 7riVA vertices, for which we choose a dipole form, 

u{k) = WTW ■ (17) 

The numerical calculations are performed with a characteristic momentum cut-off scale 



A = 0.8 GeV, just a little softer than the measured axial form factor of the nucleon |52|, |53 
- although the results are relatively insensitive to the precise value of A, as illustrated 
below. The ratio of the 7rA^A to ttNN couplings can be determined from SU(6) symmetry 
(g-wNA/ 'g-wNN = a/72/25), however, in the numerical calculations we consider a range of values 
for the ratio. SU(6) symmetry is also used to relate matrix elements of the twist-2 operators 
in the bare A and A^-A transition to those in the bare nucleon. Lattice calculations of A or 
A^-A transition matrix elements will in future test the reliability of this approximation. 



2 While the heavy baryon limit applies strictly when << M, the form factor, u(k), strongly suppresses 
all of these integrals for m, above 400 MeV and thus the heavy baryon expression provides an adequate 
description of the meson loops in the region where they are large and rapidly varying. 
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FIG. 3: Contributions to the wave function and vertex renormalization of the nucleon matrix 
elements of the operators Of 1 '" i = q,Aq,5q, in Eq. (g). Solid, double and dashed lines denote 
nucleon, A and pion propagators and the crossed circle and box indicate the insertion of the relevant 
operators. Diagrams and denote the contributions to wave function renormalization (a 
derivative with respect to the external momentum is implied). 

The renormalization constants for the spin- independent, helicity and transversity opera- 
tors are given by 

Z q l = l + Z™ + Z$ + Z\$, (18a) 
ZaI = 1 + Z™ + Zj* + Z$ + Z*$ + Z£* + Z™ T + Z*J* , (18b) 

y-l 1 1 yNN 1 yNA , yAN , 7 AA , ^ /vtad , ^ yNWT , ^ 7AWT /1 o \ 

Z <5q = 1 + ^1,P + A,P + A,P + Z 1,P + 2 1> P + 2 1,P 2 ' P ' 

The contributions from the coupling to nucleon intermediate states are given by: 

yNN _ 9a f°° k 4 u 2 (k)dk 

and 

yNN 1 9 2 A r k'u\k)dk 



3 (4tt^)2 7 ^(fc) 

for the unpolarized and polarized operators, respectively. One can explicitly verify that the 
LNA behavior of these contributions is m^logm 2 . The A contributions to the unpolarized 
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and polarized operators are equivalent, 

yAA _ yAA __ ^0 9nNA 9 A f°° k 4 u 2 (k)dk 

W l ' P " 21gl NN {^Uf] u(k)(u(k) + AM) 2 1 ) 

while the NA transition contributes only to the spin-dependent operators, 

yAN _ yNA ^6 9ttNA 9 a f°° k 4 u 2 (k)dk , . 

These contributions are illustrated in the middle row in Fig. || Expanding these terms for 
small m n , one finds that the leading non-analytic terms associated with the A and N-A 
transition contributions enter at orders m 4 logm 2 and m^, respectively. The contributions 
from the tadpole diagrams, which are independent of gA, are also identical for the unpolarized 
and polarized cases, and given by 

Had _. ^tad _. 2 f°° k 3 u 2 (k)dk 



The tadpole contributions also enter at order m 2 \ogrn 2 . ||43|| , as can be verified directly from 
Eq. (H). 

While the inclusion of the A resonance is important for quantitative descriptions of 
baryon structure, we also know from phenomenological studies that the higher order (in 
m n ) Weinberg- Tomozawa contact term |54|, [55| plays a vital role in low energy S'-wave pion- 



nucleon scattering ||56|| . Because of the Adler-Weisberger relation |57j between nN cross 
sections and g^, any term which affects ttN cross sections should also have some effect on 
9a ||58|| . In fact, within the CBM Morgan et al. [|58| found that this term largely resolves 
the discrepancy between the bag model value of gA = 1-09 and the empirical value of gA 
for bag radii R G (0.9, 1.1) fm. In the present treatment, since we do not use the CBM 
explicitly, but rather parameterize the pion source via the phenomenological form factor 
u(k), we determine the overall strength of the Weinberg- Tomozawa term so as to reproduce 



the contribution found in the CBM, as outlined in Ref . |58[ . The relative contributions of the 
diagrams with iV and A intermediate states, however, illustrated in the last row of Fig. || 
can be fixed by SU(6) symmetry. These contributions to the operator renormalization can 
then be written: 

7 ivWT n f°° k 4 u 2 (k)dk 







u 2 (k) 



7 awt _ Cwt n rx\ 9nNA f°° k 4 u 2 (k)dk 

Zl > P "^ (+V) ^7o W (fc)Mk)+AM) ' (25) 
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for the N and A intermediate states, respectively, where Cwt is the overall normalization. 
For the above range of R, the physical value of g& can be reproduced to within a few 
percent for the corresponding range of Cwt €= (0.21,0.30). In the following numerical 
analysis, we use this range as an estimate of the systematic error on the Weinberg- Tomozawa 
contribution. Even though the non-analytic behavior of the integrals in Eqs. ( ]2jj ) and ( 25 ) 
is m% or higher, their contributions are found to be significant. Note, however, that the 
Weinberg- Tomozawa terms contribute only to spin- dependent matrix elements, and make 
no contribution to unpolarized matrix elements. 

With the exception of the matrix elements of the unpolarized, n = operator, the renor- 
malization of each moment of the various distributions is independent of n. The n = 
operator, which corresponds to the isospin charge, is not renormalized — additional contri- 
butions from operator insertions on the pion propagator cancel those shown in Fig. |3|. 

The pion mass dependence of the various contributions to the wave function and oper- 
ator renormalization is shown in Fig. |j. For the ratio of the couplings, g^NA/ gnNN, SU(6) 
symmetry is assumed. The relative size of the terms Zf N and Z AA in the spin-dependent 
and spin-independent cases already makes it clear that intermediate states involving A reso- 
nances are much more significant in the former case. In particular, whereas Z^p does little 
to counter the effect of the wave function renormalization, the A contributions Z A p, Zfp, 
and Z^p essentially cancel its effect. 

To explore the sensitivity of the results to the strength of the A contribution, in Fig. |5| 
we show the combined effect of the pion dressing on spin-averaged (upper panel) and spin- 
dependent (lower panel) nucleon matrix elements in Eq. (|U|) for a range of values of the 
ratio g n NA/ 9-kNN- For illustration we choose values of g w NA/ 9wNN equal to zero (no A 
states), ^72/25 (SU(6) coupling) and 2 (phenomenological value needed to reproduce the 
width of the physical A resonance). In the unpolarized case, the effect of this variation is 
relatively small — less than 3% over the entire range of pion masses considered here. In 
contrast, the effect of the A on the helicity and transversity moments (matrix elements of 
the spin-dependent operators Q^"-^ n an d O^'"^™) is far more significant. If the contribu- 
tion from the A (and the N-A transition matrix elements) is ignored {g n NA/ 9-kNN — 0), 
ZiIZaA phys = 0.90, while including these contributions with SU(6) couplings increases 
this to Z 2 /Za | phy S = 1.01, and to Z 2 /Z/± g \ P h ys = 1.04 at the phenomenological value 
g-KNA/ 'g-wNN = 2. Consequently, when the effects of the A are included with a coupling 
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FIG. 4: Contributions to the pion loop renormalization of the matrix elements of the twist-2 
operators required to evaluate the moments of the PDFs. The upper panel shows nucleon wave 
function renormalizations (Z^ , Z%) and spin-independent operator renormalizations. The lower 
panel shows the contributions to the renormalization of spin-dependent operators, and the shaded 
region is an estimate of the uncertainty in the Weinberg-Tomozawa term, Z^p = Zfp^ T + Z^p^ T . 
The flv NA I 9-k N N coupling constant ratio is set to the SU(6) symmetric value of yw2/25. 
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FIG. 5: Pion dressing of the matrix elements of the spin-independent (upper panel) and spin- 
dependent (lower panel) operators in Eq. (|3|) for various values of the ratio of coupling constants, 
9ttNa/9ttNN- The shading in the lower panel indicates the variation about the dashed curve 



(9ttNa/9itNN = y/72/25) caused by the uncertainty in the Weinberg-Tomozawa term. The be- 
havior of Z^jZ^n is similar to that of Z2/Z& q . 
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FIG. 6: Dependence of the renormalization of the nucleon matrix elements in Eq. ( |l3| ) on the 
dipole mass parameter, A, at the physical pion mass, for the SU(6) value of gnNA/gnNN- 

constant which is consistent with phenomenology, one finds that there is almost no curva- 
ture in the extrapolation of the spin-dependent moments. This result is relatively stable 
against variations [|5l| |53| in the dipole mass parameter, A, in the range ~ 0.7 — 1.0 GeV - 
especially for the spin-dependent moments, as illustrated in Fig. |6|. 

Matrix elements of the twist-2 operators (|3|) in bare nucleon states will necessarily be 
analytic functions of the quark mass (m q ~ m^), so the one pion loop renormalization 
described above is the only possible source of LNA contributions. Consequently, the LNA 
behavior of the matrix elements in Eq. ( |13| ) will be given by 

(<T- M ")lna = Zf A - Z^ A , i = q, Aq, Sq . (26) 

If the N-A mass splitting is artificially reduced to zero, A intermediate states become de- 
generate with the corresponding nucleon intermediate states, and the respective A diagrams 
formally give rise to LNA contributions. Leaving the g n NA/ 9-kNN ratio free, the coefficients 
of the LNA contributions (the log term) to the various matrix elements can then be 
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written: 



W" M ">LNA = 
\y 5q /LNA 



^•LNA 



^LNA 



7 LNA 7 LNA 
^2 ~ ^Aq 



1 



3 + 



7LNA 



7 LNA 
^Sq 



1 



2 - 
2 - - 



16 gjwVA 

27 
4^ 



NN 



NA 



9 9ttnn 

4 gl N A 

9 SttNN 



«A+1 



(27a) 
(27b) 
(27c) 



(4vr/.) 2 

This makes it clear that, whereas an increase in g n NA/gnNN from (no A contributions) 
tends to increase the effective coefficient of the chiral logarithm in the unpolarized case, for 
the spin- dependent operators it acts to suppress it. Indeed, assuming the bare axial coupling 
qa = 1.26, at g-rtNA/ 9ttNN = 2.43 the LNA coefficient for the polarized moments is zero, and 
for larger values it even becomes positive. Whilst this exact cancellation is an artifact of 
setting AM = 0, it highlights the significant role played by the A resonance. 

From this analysis and the numerical results shown earlier, one can conclude that the 
inclusion of the A resonance will cause only a minor quantitative change in the extrapola- 
tion of unpolarized moments, and in practical extrapolations of lattice data the A can be 
neglected with no loss of accuracy, given the current uncertainties in the data. In contrast, 
the A leads to a qualitatively different picture for the extrapolation of the spin-dependent 
moments and must be included there. 

There are a number of possible approaches that can be taken to account for the A 
contributions. One strategy would be to include the one-loop renormalizations numerically 
in the extrapolations, along the lines of the calculation of self-energies in the hadron mass 



extrapolations in Refs. |39|, One could also replace the momentum integrals in the 

expressions for Z\ and Z 2 with discrete sums over momenta which are available on the 
lattice, J d 3 k — > (l/V)(27r/a) 3 J2k x k y k z i wri ere V is the spatial volume of the lattice, as in 
the analysis of the p meson mass in Ref. |S9j (see also [§t]]). Because of the discretization 
of space-time on the lattice, the lattice momenta are restricted to values = 27Tn fl /aL fJi , 
where is the number of lattice sites in the /i direction and the integer runs between 
— and +L fl /2. We have checked that at large m 2 the differences between the integral 
and discrete sum are only a few percent or less, however, at small m 2 values the momentum 
gap between = and the minimum momentum allowed, k^ = ±2n/aL fJ- , may introduce 
corrections. 

Although this procedure is more accurate in principle, in practical extrapolations of lattice 
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data it is not as straightforward to implement as an extrapolation formula based on a simple 
functional form would be. For this purpose it is more useful to preserve the simplicity of a 
single formula which interpolates between the distinct realms of chiral perturbation theory 
and contemporary lattice simulations, as proposed in Refs. [O, [IB] . In order to test whether 
one can continue to apply a modified form of the extrapolation formula in Eq. (|5J) to lattice 
data for the spin-dependent moments, as well as the spin- independent, we attempt to fit the 
pion mass dependence of the renormalization factors shown in Fig. ||] using the form 

,2 



ry I ry . n 2 i 2 1 _ 

z 2 /Zi = on + + Ta t \? m 7r lo s 



ml 



q,Aq,5q } (28) 



(4tt/ w ) 2 ' 

with «j, Pi and /ij treated as free parameters, but with ^lna fixed to the values obtained 
analytically in the limit AM — > 0, as shown in Table [II]. The fits to Z 2 /Zi (i = q, Aq, Sq) 
are illustrated in Fig. [7| for the average of the g n NA/ g^NN values from SU(6) symmetry 
(yj 72/25) and phenomenology (2). Fits for other values of the coupling are equally good. 



It is remarkable that the LNA form (^) is indeed able to reproduce the full calculations of 
Z 2 /Zi with such high accuracy, given that the full calculations include higher order effects 
(in m^) associated with the A and Weinberg- Tomozawa contributions. The best fit values of 
/i, shown in Table [Tl], are only slightly smaller than those found in earlier work M. Note 
that the functional form in Eq. (^) does not include the modifications designed to ensure 
the correct heavy quark limit, as in Eqs. (||H(fJ) — incorporating this constraint leads to 
only marginal changes in the parameter [i |T6| . 

As discussed above, excluding the isospin charge, all moments of each operator are renor- 
malized in the same manner. Hence, our conclusions regarding the inclusion of the A-isobar 
apply equally well to extrapolations of qa = (l)Aw-Ad and all other moments of the helicity 
and transversity distributions. 

IV. EXTRAPOLATION OF LATTICE DATA 

Having established that the LNA formula, Eq. (|28|), provides a good approximation to 
the full calculation, in this section we examine the effects of extrapolating the available 
lattice data on the twist-2 PDF moments using the forms @, (§) and (H), with the LNA 
coefficients determined in the limit AM — > 0. Rather than show the moments versus the 
scale and renormalization scheme dependent quark mass, Figs. || || and |10] give the moments 
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0.32 
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25 #A + 1 


0.39 




0.25 


i8„2 , 1 
25#A + 2 


0.29 


1.85 


A.5lg 2 A + 1 


0.38 


0.48<& + 1 


0.25 


0.48 5 i + i 


0.30 


2 


T#A + 1 


0.37 




0.24 


2 n 2 4- 1 

g^A + 2 


0.29 



TABLE II: Fits to the dependence on m 2 of the calculated renormalization factors, obtained by 
varying ctj, 0i and in Eq. (|2~8l). The LNA coefficients and the mass parameters /ij are shown for 
various values of the g n NA/ '9-kNN ratio: (no A), \/ 72/25 (SU(6)), 2 (phenomenological value) 
and 1.85 (average of SU(6) and phenomenological values). 
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FIG. 7: Fits to the calculated renormalization factors Z2/Z1, i = q,Aq,5q, in Eq. (|T^) using the 
functional form in Eq. (p8|), as a function of m 2 . The ttNA to nNN coupling ratio has been set 
to the average of the SU(6) and phenomenological values (g n NA/gnNN = 1.85), and gA to the tree 
level value, 1.26. The corresponding values of /i are given in the third row of Table O. 
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of the unpolarized, helicity and transversity distributions, respectively, as a function of the 
pion mass squared. The data have been extrapolated using a naive linear extrapolation 
(short-dashed lines), as well as the improved chiral extrapolations with the LNA chiral 
coefficients and values of \x given in Table O, with m& )Tl fixed at 5 GeV [03] . 

For the spin-dependent moments, four curves are shown in Figs. |9| and [OJ: the long- 
dashed curves correspond to ignoring A intermediate states (g^NA = 0), while the central 
solid lines in each panel of the figures include the effects of the A with a coupling ratio 
9wNa/ QttNN = 1-85 (average of SU(6) and phenomenological values) and the central value 
for the Weinberg- Tomozawa coefficient, Cwt = 0.255. The upper and lower solid lines 
correspond to (JttNa/ 9nNN = 2, C WT = 0.30 and q^na/ 9ttNN = a/ 72 / 25 ' c wt = 0.21, 
respectively. Because the effect of the A is almost negligible for the unpolarized moments, 
these curves are all essentially collinear, and for clarity only one is shown in Fig. || The 
extrapolated values are shown in Table |T|, along with the associated errors (which are 
described in the Appendix) and the experimental values for the unpolarized and helicity 
moments |4], |^, ^|, [7|] (there are currently no data for the transversity moments). Note that 
even though there is no curvature expected for large m n , the slopes of the linear and LNA fits 
differ slightly at large m n values due to the constraints of the heavy quark limit incorporated 
into the forms ([3J), ([|) and (|]). 

With respect to the moments of the unpolarized PDFs, this analysis confirms our earlier 
finding that it is essential to incorporate the correct non-analytic behavior into the chiral ex- 
trapolation. When this is done, there is good agreement between the extrapolated moments 
at the physical pion mass and the corresponding experimental data. On the other hand, for 
the polarized PDFs we have the surprising result that once the A resonance is included, the 
effect of the non-analytic behavior is strongly suppressed, and a naive linear extrapolation 
of the moments provides quite a good approximation to the more accurate form. 

In the case of the axial charge (the n = moment of Am — Ad), the extrapolated value 
lies some 10% below the experimental value, with an error of around 5%. However, 
appears to be particularly sensitive to finite volume corrections, with larger lattices tending 
to give larger values of Furthermore, there is some sensitivity to the choice of action 
- simulations with domain wall fermions (DWF), which satisfy exact chiral symmetry, 
are found to give larger values than those with Wilson fermions [^BJ. As almost all of 
the currently available lattice data are obtained from very small lattices (L ~ 1.6 fm), we 
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FIG. 8: The lowest three non-trivial moments of the unpolarized distribution u — d, extrapolated 
using a naive linear fit (dashed lines) and the improved chiral extrapolation (solid lines). The stars 
indicate the experimentally measured moments at the physical pion mass, and the lattice data are 
taken from the sources listed in Table |, where the various plotting symbols are defined. 



26 



1.60 



1.40 



< 

k 1-20 

A 

v 1.00 

n 

^0.80 



0.60 




0.2 



0.40 



^ 0.30 

i 

A 

v 0.20 



0.70 
0.0 



0.30 



0.2 



0.20 



i 

< 

A 



? 0-^0 



0.00 



0.0 0.2 



0.4 



0.6 



0.8 




0.4 



0.6 



0.8 



0.4 0.6 0.8 



1.0 



1.0 



1.0 



• — ♦ 



1.2 



1.2 







i j 


--#--1-- \ * — 

I T T T i 



1.2 



FIG. 9: The lowest three moments of the helicity distribution Au — Ad, extrapolated using a naive 
linear extrapolation (short-dashed lines) and the improved chiral extrapolation described in the 
text. In each panel, the long-dashed lines correspond to fits with no A and the LNA coefficient 
determined from xPT, while the solid lines are fits obtained using g^Nk/ 9irNN = 2 (upper solid 



curves) and y72/25 (lower solid curves). The lattice data are taken from the sources listed in 
Table |. 
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Moment Value Extrapolation errors a n , Aa n ,6a. 





Experimental 


Extrapolated 


Statistical 


A, WT states 






{x) u ~d 


0.145(4) 


0.176 


0.012 


0.0008 


0.022 


0.141 


(X 2 ) u -d 


0.054(1) 


0.054 


0.015 


0.0003 


0.007 


0.044 


{X 3 )u-d 


0.022(1) 


0.024 


0.008 


0.0001 


0.003 


0.019 


WAu~Ad 


1.267(4) 


1.124 


0.045 


0.020 


0.022 


1.084 


(x)Au-Ad 


0.210(25) 


0.273 


0.015 


0.005 


0.005 


0.262 


(x 2 )Au-Ad 


0.070(11) 


0.140 


0.035 


0.003 


0.003 


0.135 


(l)8u-6d 




1.224 


0.057 


0.019 


0.025 


1.187 


(x)su-8d 




0.506 


0.089 


0.008 


0.010 


0.490 



TABLE III: Values of the unpolarized, helicity and transversity moments, extrapolated to the 
physical pion mass using Eqs. (||), (§) and @ and the parameters in Table [0]. The experimental and 
systematic errors listed here are described in the Appendix. For comparison, experimental values 
of the moments where known (unpolarised values from moments of distributions of Refs. 0, ^, 
helicity moments from Ref. |?J in scenario I (NLO)) and the best fit parameters (a n , Aa n , 5a n ) are 
also listed. 

consider the current level of agreement quite satisfactory. 

Additionally, there is some uncertainty arising from the inclusion of the heavy quark limit 
in our fits; if this constraint is omitted, the large m 2 behavior of our fits coincides with the 
linear fits that are shown as one would expect. For (x)s u -Sd ((x) Au-Ad), a fit ignoring the 
heavy quark limit gives a physical value of 0.559 (0.257) rather than 0.506 (0.273) as given 
in Table HJ (with a smaller effect in the other moments). 

The uncertainty in the experimental determination of the higher moments of the spin- 
dependent PDFs is considerably larger, and from the current data one would have to conclude 
from Fig. ^ that the level of agreement between experiment and the extrapolated moments 
is acceptable. Clearly the scatter in the lattice data for the second moment means that at 
present we cannot have much confidence in the predicted value. We do note, in addition, 
that this is one case where there is a tendency for the full QCD points to lie somewhat below 
the quenched QCD results. It is obviously of some importance that this issue be resolved in 
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FIG. 10: The lowest two moments of the transversity distribution bu — 5d. All curves are as 
described in Fig. ||[ 

future lattice simulations. 
V. CONCLUSION 

The insights into non-perturbative hadron structure offered by the study of parton dis- 
tribution functions makes this an extremely interesting research challenge. It is made even 
more important and timely by the tremendous new experimental possibilities opened by 
facilities such as HERMES, COMPASS, RHIC-Spin and Jefferson Lab. Lattice QCD offers 
the only practical method to actually calculate hadron properties within non-perturbative 
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QCD, and it is therefore vital to test how well it describes existing data. Because current 
limitations on computer speed restrict lattice simulations to quark masses that are roughly a 
factor of 6 too large, one must be able to reliably extrapolate the lattice data to the physical 
quark (or pion) mass in order to compare with experiment. 

Traditionally such extrapolations have been made using a naive linear extrapolation as a 
function of (or quark mass). In Ref. [TE\, Detmold et al. showed that it was essential to 
include the leading non-analytic behavior of chiral perturbation theory in this extrapolation 
procedure. Only then were the existing lattice data for the moments of the unpolarized 
parton distribution functions in agreement with the experimental moments. Here we have 
confirmed this conclusion by calculating the next-to-leading non-analytic behavior within a 
chiral quark model, including the A-isobar, and showing that it led to precisely the same 
conclusion. 

We have also investigated the variation of the moments of the polarized parton distribu- 
tions to next-to-leading order. In this case the inclusion of the A-isobar makes a dramatic 
difference. Indeed, once the A is included, the helicity and transversity moments show little 
or no curvature as the chiral limit is approached and a naive linear extrapolation formula is 
reasonably accurate. In case a more accurate extrapolation procedure is desired, we propose 
convenient formulae which suitably build in the non-analytic behavior in both the unpolar- 
ized and the polarized cases. The value of Qa extracted from the extrapolation procedure 
at the physical pion mass is within 10% of the experimental value. Given the sensitivity 
of this quantity to lattice volume (current simulations use quite small lattices) and quark 
action (domain wall fermions tend to give a larger value of qa than Wilson fermions), we 
consider this a very satisfactory result. We look forward with great anticipation to the next 
generation of lattice simulations of parton distribution functions at smaller quark masses 
and on larger volumes. 
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APPENDIX A: STATISTICAL AND SYSTEMATIC ERRORS 

In this Appendix, we describe the estimates of the statistical and systematic errors in our 



fits that are presented in Table [11 1 



To determine an estimate of the error associated with the statistical uncertainty of the 
lattice data, we use the estimated standard deviation. For data, fa, and weights, cjj, given at 
abscissae X{ (i — 1, . . . , n), the estimated standard deviation of a fitting form f(x; a) with 
parameters a is: 



0o 



1 n 



Xi;a )) 2 , (Al) 



where c?o are the best fit parameters. The statistical errors assigned to the fits are then 
determined by varying the fit parameters (a n , Aa n , 5a n ) from their optimal values (given in 
the right-most column of Table |T|) to obtain an increase of unity in the standard deviation. 

In order to estimate the systematic errors arising from the form of our fits, we first 
consider the uncertainty in the values of g^Nk/ g-nNN an d Cwt, taking half the difference 
between the physical values of the moments obtained with g^NA/ 'g-wNN = 2, Cwt = 0.30 and 
g-K TV a / giv n n = a/72/25, Cwt = 0.21. We also consider the uncertainty in the fit parameter 
\i by taking half the difference between the physical moments obtained with \i 20% above 
and below the fits obtained in Table 0. The resulting systematic uncertainties are listed in 
Table IJT]. 
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